On the running of the bare coupUng in SU(N) 
lattice gauge theories 



Chris AUton", Michael Teper^ and Aurora Trivini" 

''Department of Physics, University of Swansea 

^ Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford 0X1 3NP, U.K. 



Abstract 

Interpreting the way that the SU(3) bare lattice coupling runs with the lattice spacing is 
complicated by the fact that there is a smooth cross-over region in which the strong coupling 
expansion transforms into a weak-coupling one. For > 5, however, there is a first order 
bulk transition that cleanly separates the strong and weak coupling regimes. We find that 
in this case the calculated string tension can be readily fitted throughout the weak coupling 
region by a standard 3-loop expression modified by lattice spacing corrections of the expected 
form. While our fits demand the presence of the latter, they do not constrain the perturbative 
coupling scheme enough to enable us to extract a usefully accurate value of a(/3) in units 
of Ajjg. To resolve this ambiguity we turn to SU(3) where we use the Schrodinger Func- 
tional coupling scheme to extract a value of tqAsf as a benchmark. We then find that the 
Parisi mean-field improved coupling scheme closely reproduces this result. We also develop a 
comparison between different schemes that does not rely on the calculation of any physical 
quantity and which can therefore be applied much further into weak coupling. Again the 
Parisi scheme is favoured over the others that we compare. Using the mean- field scheme we 
have fitted the values of the string tension a^a that have been calculated for 2 < A^ < 8, to 
obtain Ajts/V^^ = 0.503(2) (40) + 0.33(3) (3)/A^2 AT > 3, where the first error is statistical 
and the second is our estimate of the systematic error from all sources. 



1 Introduction 

Consider SU(A^) gauge theories discretised onto a hypercubic lattice of spacing a, on a 4-torus, 
with SU(iV) matrices, Ui, assigned to the hnks /, and with the standard Wilson plaquette 
action. The partition function is 



where Up is the ordered product of the SU(iV) matrices around the boundary of the plaquette 
p. (Although we shall be using the plaquette action in this paper, a parallel analysis could 
be carried out for any other lattice action.) The parameter f3 in the lattice partition function 
is proportional to the inverse lattice bare coupling. This defines a running coupling on the 
length scale a, in what is often called the lattice scheme: 



In practice one frequently wishes to compare physical quantities that have been calculated 
over a similar range of /?, but not at precisely the same values. For example the Sommer 
parameter PQ or the deconfining temperature Tc [2] in terms of the confining string tension 
a. One then needs to interpolate such quantities to common values of f3. It is well-known 
that simple perturbative interpolations will not fit the calculated values. Thus one is typically 
reduced to using, for example, interpolating polynomials [3] that bear no relationship to the 
weak-coupling expansion, or, if one tries to use a power series in g1 (x l/f3, one finds that the 
higher order fitted terms have large coefficients of oscillating signs, so that if one extrapolates 
further into weak coupling, the values diverge away from the expected 3 loop form until truly 
asymptotic values of the scale 1/a [1]. This is, of course, a symptom of the well-known fact 
that in the range where one currently performs calculations, no plausibly simple perturbative 
expression for g1{a) even begins to be adequate. This is unfortunate because it means that one 
cannot exploit dimensional transmutation in a simple way, to transform the value of the lattice 
bare coupling to a statement of how large the lattice spacing a is in units of the corresponding 
A parameter, which can then be expressed in terms of, say, Ajjg in the theoretically and 
experimentally well-studied MS coupling scheme. (See [5] and |6] for recent reviews.) 

This old problem has been approached both by stressing the need for improved lattice 
coupling schemes [3, [8] and the need for lattice spacing corrections [9j. Determining what 
is important is, however, rendered ambiguous by the fact that there is a smooth cross-over 
between strong and weak-coupling in SU(3), so one does not really know from which value of 
(3 it is appropriate to attempt a weak-coupling expansion in powers of and a. 

In this paper we use the fact that for SU(A^ > 5) there is a first order 'bulk' transition ^10] 
separating the weak and strong coupling ranges, to remove the ambiguity of where one might 
expect a weak coupling expansion to be applicable. This enables us to quantify the importance 
of retaining O(a^) lattice corrections in addition to the usual continuum perturbative variation. 





2N 



(2) 




1 



While our SU(6) and SU(8) calculations prove accurate enough to establish the need for 
O(a^) corrections in the relationship between a and g'^{a), the range in a is not large enough 
to usefully constrain the coupling scheme and hence the value of aAj-j^. In fact experimental 
studies are not accurate enough either. (See for example Fig.lO in [llj.) Fortunately there 
exists in SU(3) an accurate calculation of the running coupling in the 'Schrodinger functional' 
(SF) scheme that covers an energy range comparable to that of experiment, i.e. up to ~ Mz, 
and with appreciably smaller errors [121 [E] • We shall use this scheme to obtain from the values 
of a/ro calculated in [3l [H] the continuum value of tqAsf and hence of roAjj^. We compare 
this to what one obtains with some improved coupling extrapolations, and find that the Parisi 
mean-field improved coupling scheme closely matches the SF result. We simultaneously 
perform a comparison with the SF scheme that does not involve the calculation of any physical 
quantity and therefore can be carried out to much weaker coupling. This also points to the 
'goodness' of the mean-field scheme. Motivated by this we use this scheme for iV 7^ 3 to obtain 
continuum values for A-^jg/ ^/a for all A^, and in particular for — 00. 

In the next Section we discuss the weak coupling behaviour of a on g'^{a) in more detail, 
and establish our notation. We then move on to our calculations and then attempt to quantify 
the systematic errors on our final results. We end with a brief summary and our conclusions. 

An abbreviated version of this work has been presented at Lattice 2007 [16]. This study 
forms part of a more detailed and extensive analysis that will appear elsewhere [17] . 



2 The lattice running coupling 

Ideally we would wish to be able to determine from the value of g\{a) = 2N/(3 what is the 
lattice spacing a as expressed in units of the corresponding physical scale Al (which can 
then be converted into a more familiar scale, such as A;^, using a one- loop calculation [18j). 
Although we know this to be possible in principle (dimensional transmutation) in practice 
doing so for the SU(3) gauge theory (and QCD) has proved notoriously ambiguous, despite 
the fact that we know the /3-function in the lattice scheme to 3 loops [19]. There are some 
plausible reasons for this which will become apparent in the following discussion. 

To establish how g1{a) is related to a we determine how a varies with g\{a)- To do this we 
express a in units of some physical mass ix that we calculate in lattice units, i.e. as afi{a). Now 
we know that different choices for the quantity ^ will have different lattice spacing corrections: 

where we have chosen to use the string tension a to set the scale of these corrections. As 
an example, the ratio mQ++/y/a, where mo++ is the hghtest scalar glueball mass, has Co- ~ 2 
to 3 depending on A^, while Tc/y/a, where Tc is the deconfining temperature, has ~ 1/3. 
Thus it is clear that the variation of afi{a) with g1{a) must include a dependence on a that 
has the functional form shown in eqn([3]), in addition to the variation expected for a running 
coupling in the continuum theory. Generically we should expect Co- = 0(1) as illustrated by 
the examples just quoted. There is a similar dependence that follows from the fact that the 
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/?-function for g1{a) will have lattice spacing corrections 

dgl 



d log 



f3o9l + f3i9i + f3^9l + ... + 0{a'). (4) 



Here the Pi are the coefficients of the continuum /3-function; Pq and Pi are scheme-independent 
while P2 requires a 2-loop calculation relating gj^ to a coupling for which P2 is already known 
[T9] . Thus, as has been emphasised in [9J, the dependence of a on (/I, (a) will need to incorporate 
such lattice spacing corrections in addition to the perturbative expression that one obtains in 
the continuum. Using the string tension a for our scale this leads to 

, ^-4.2+(|f-||)rf«(,S) 

where we make explicit the quantities that depend on the coupling scheme s being used (here 
s = L). We use the standard definition of the A parameter, where the constant term in the 
expansion of l/g'^{a) is absorbed into the scale of the logarithm. In eqn([5]) the terms that 
involve only Pq and Pi constitute the exact 2-loop continuum result. (That is to say, it is the 
exact result when Pj>2 = 0.) For higher orders we are not aware of any such neat closed form 
in terms of elementary functions, so we present their contribution as a power series in g^. Note 
also that although the coefficient is a power series in g^, we shall, following usual practice, 
treat it as a constant in our ffis, since g1 does not vary very much in the region where the 0{o?) 
correction is significant. Finally we remark that we could introduce the lattice corrections in 
different ways. A small change would be to use, say, the mass gap ma rather than the string 
tension a to set the scale for the 0{o?) corrections in eqn([n]). A more radical change would 
be to substitute for each occurence of a^a on the RHS of eqn([5]) the expression provided by 
eqn([5]) itself. After iteration this would transform the lattice spacing corrections from a power 
series in a^a into a power series in the perturbative factor in eqn([5]). All such variations are 
a priori equally valid. 

In addition to these lattice spacing corrections, a further complication is that the 5'l(ci) 
coupling scheme is expected to have large higher order perturbative corrections. This follows 
from the large ratio between A^ and I^-ms [18] : 



This implies that in the relationship between the two couplings, g'^-^ = g\{l + 'jg'l + 0{g'l)), 
the coefficient 7 must also be large, since one can easily see that Aj^/Al = exp(7/2/3o). So if 
g^^ is a 'good' scheme with modest higher order terms in the /3-function, which is something 
we shall assume from now on (see [20]), this will almost certainly not be the case for g1. This 
is con&med by explicit calculation of the three-loop coefficient [19] where one finds 

P^ 6.9/3p (7) 
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Since we only know the /^-function to 3 loops, it would be wise to seek a lattice coupling 
scheme where there is less reason to expect large higher order corrections. This problem has 
of course been appreciated for a long time and there have been extensive efforts to improve 
the lattice coupling (for a review see [8].) Perhaps the simplest and oldest suggestion is the 
'mean-field' improved coupling of Parisi [7] 

which has a nice physical motivation as the effective coupling experienced by a background 
field (in a simple approximation) [7j. Since the expansion of (Trf/p) in g1 is known to 3-loops 
pi] we can substitute 

1 = \{^TtU,) = \{l- u,gl - u,gl - u,gl + 0{gl)) (9) 
9i 9l 9l 

into the 3- loop /?-function for g1 to obtain the 3-loop /^-function for g], giving 

P'2 = P2 + - A^i "^^^ 2.3/3p (10) 
which indeed promises smaller higher-order corrections than eqn([7]). Similarly one finds 

Am ^ A^A^ ^ ^ 38.853exp (-^1 ^ 2.633 (11) 

using eqn(l6]) and ui = {N'^ — 1)/8A^, which again indicates that the mean-field improved 
coupling scheme is a 'good' one. Motivated by this we shall focus on this coupling scheme, 
and some variations thereof, in our later calculations. 

Although the need to improve the lattice coupling has long been recognised, the suggestion 
[9] that O(a^) lattice corrections are also needed has in practice been more controversial. As 
demonstrated above, there can be no question about the presence of such corrections. But one 
might question their importance. How plausible is this? As remarked above, we expect that 
generically = 0(1) in eqnQ, and since a^/o^{a) and other masses are typically calculated 
to a precision ^ 1%, we would need a^cr ^ 0.01 for the O(a^) correction to be negligible. In 
SU(3), for example, a^a ^ 0.01 at /3 = 6.5. This coupling region is at the upper edge of the 
range in which useful calculations are usually performed. Thus generically we should expect 
O(a^) corrections to be important in current calculations. Of course it might have been that 
the coefficient just happened to be unexpectedly small, but the calculations in this paper 
will show that this is not the case. 

Eqn([5]) is only valid in weak coupling. This is as true of the lattice corrections as of the 
perturbative factor. In strong coupling the appropriate relationship is a very different one. 
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(See ^22j for a detailed review.) A problem that now arises for small A^, and in particular for 
SU(3), is that there is a smooth cross-over between the strong and weak coupling regimes, 
with a mid-point that occurs close to the value of (5 where one typically starts to calculate 
interesting physical quantities. This is illustrated in Fig. [1] where we show how ay/ a varies 
with /? for SU(3). The point of inflection, around /? G [5.50, 5.60], provides an estimate of the 
mid-point of the cross-over. Clearly there must be some range of (3 beyond that where the 
functional form is still not pure weak coupling, but a priori we do not know how far that range 
extends. This complication disappears for > 5 where the strong-weak coupling transition 
becomes a first order phase transition. This is illustrated in Fig. [2] for SU(8). Here we see a 
large hysteresis indicating a strong first order transition. The simultaneous large jump in the 
average plaquette, and hence action, shows this to be a 'bulk' transition where the physics 
changes on all length scales. It is reasonable to expect that on the weak-coupling side of this 
bulk transition there are no strong-coupling artifacts. We shall take advantage of this feature 
by first testing our weak-coupling fit in eqn([S]) to the values of a^a obtained on the whole 
weak coupling branch of the SU(8) theory, and then decreasing A^ to A^ = 3 (and N = 2) 
having established in the much less ambiguous large- A^ case what kind of fit is needed. 

Our assumption that there are no strong-coupling artifacts on the weak-coupling side of the 
'bulk' transition is not only plausible but is supported by analytic calculations in a related but 
simpler context: SU(A^) lattice gauge theories in 1 + 1 dimensions. Here (with the plaquette 
action) there is a strong-weak coupling cross-over at finite A^ that becomes a third-order phase 
transition at A^ = cxo [30]. In this limit one can calculate the string tension as a function of 
j3 analytically [30] and one finds a simple expansion in powers of with no strong-coupling 
artifacts anywhere on the weak-coupling branch. While this does not prove that the same will 
be true in the case of = 3 + 1, it provides an additional argument that such artifacts will 
be either absent or highly suppressed. 

An interesting question concerns the functional form, for A^ < 4, in the extended cross-over 
region where the strong coupling expansion gradually transforms into a weak-coupling one. 
We can obtain some intuition from the D = 1 + 1 case where the whole problem reduces to a 
single SU(A^) integral [30]. In particular, for SU(2) the running of the coupling becomes quite 
elementary, 

J^^dzz{l — z'^)^ exp {(3z) 

Up = exp{— a cr} = — J ; (13) 

dz{l — ^2) 2 exp {(3z) 

and one can analyse the behaviour in the cross-over region. We note that the terms in eqn (|T3l) 
are obtained by applying d/d/3{l — or (1 — d'^/dfS'^)^ to 

1 

/ dz exp (pz) = (e^ - e"^) . (14) 
J-i P 

So on the weak coupling side we will have, in addition to the expected weak coupling expansion 
in powers of 1//5, a correction term that is an expansion in powers of e~^^, whose coefficients 
are themselves series in powers of 1/(3. It is these exponential terms that are our strong- 
coupling artifacts. They arise from the compact integration range of the plaquette in eqn flT^ . 
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and this general origin suggests that they will be present in some form in higher dimensions. 
This is most convincing for the average plaquette, where such 0(e~^^) corrections would 
create serious complications for the standard methods of extracting the gluon condensate 
in Z) = 3 + 1 since they (and their larger- homologues) would dominate over the O(a^) 
condensate contribution. We note the close connection of these corrections to ^Zn vortex- 
instantons' inD = l + l. InD = 2 + l these will generalise to Z]\f 'monopole-instantons' 
and in D = 3 + 1 to 'monopoles', as well as Z^r vortex lines and sheets. Such ultraviolet 
fluctuations will disorder small Wilson loops, and will affect the string tension in strong 
coupling. Around the cross-over one might expect some complicated contribution that is the 
ratio of powers and exponentials in jS. While it is not easy to be precise about this, it is clear 
that the functional form suggested by our D = 1 + 1 example is very different to anything one 
might imagine on the basis of weak coupling arguments, and it will be no surprise if we find 
that at small N we have to go further into weak coupling to be able to apply standard weak 
coupling fits. 

Of course, using SU(A^) gauge theories at larger to teach us something about SU(3) is 
only convincing if it is clear that SU(3) is 'close to' SU(oo). That this is so in the context of 
the bare coupling was demonstrated in Fig. 7 of |T0] using calculations for 2 < A^ < 5. Here 
we repeat the exercise with calculations that not only go to larger A^, but also to smaller 
values of a at smaller A^. So, for each value of (3 at which we calculate the string tension (see 
Section [3A1) we extract the bare mean-field improved 't Hooft coupling g]{a)N using eqn(IH]). 
We then plot it against the length scale a on which it is running, with a expressed in units of 
the string tension. Since the physics has a smooth large- A^ limit, this is a common unit up to 
corrections of 0(1/A^^) which we expect to be modest, since this is what one finds for various 
mass ratios [lOj. The resulting plot is shown in Fig. [3l This provides convincing evidence that 
the weak-coupling running of the bare coupling is in fact very similar for all values of A^. Thus 
it makes sense to extract lessons in this context from larger A^ for all SU(A^) gauge theories. 



3 The calculation 

3.1 string tensions and ro 

The string tensions used in this paper have been obtained from calculations of the ground 
state energies of confining flux loops wrapped around a spatial torus. We use the results of 
[23] supplemented in many cases by either higher statistics at the same values of f3 or new 
calculations at both higher and lower values of /3. In Tables [T]- El we list these masses, ami, 
together with the length, L, of the flux loop in each case, and the average plaquette, Up, which 
will be needed for transforming to the mean-field coupling scheme. 
To extract the string tension we use 

ami{L) = a?aL . (15) 

3/7 

The linear term has been corrected by the universal string correction [23] with a coefficient 
that corresponds to a bosonic string theory (the universality class of the Nambu-Goto string). 
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That this is in fact the appropriate universahty class for the confining flux tube in D=3+l 
non-Abelian gauge theories, has received support from many lattice calculations in both SU(2) 
and SU(3) (see for example [25]) as well as at larger [26]. The string lengths in Tables [TM] 
have been chosen to be large enough, > 3, that the string correction in eqnf llSI) provides 

only a ~ 5% shift in the calculated value of a^/o^, so higher order corrections in (the 
natural expansion parameter in the relationship between ami{L) and a^cr) should be very 
small indeed. We shall estimate the systematic error that this brings about below. 

In the case of SU(3), the calculated values of the Sommer parameter, tq, extend consider- 
ably further into weak coupling than the string tension. For this reason we will supplement 
our fits of a-y/cr with those of a/vQ, as listed in Table El These have been taken from Table 
2.6 of where, for the values at higher /3, we have translated from to using eqn(2.73) 
therein. (We note from Fig. 2. 4 of [H] that the quoted error on the conversion factor appears 
to encompass any possible dependence on a.) We remark that the lower j3 values are from |3] 
and the higher values from [15] . We have calculated the values of the action listed in Table 
These calculations have been performed on 8^ lattices and will therefore differ slightly from 
the large volume limit. However we have checked, performing a calculation on 6'^ lattices for 
the values of f3 in Table [21 that these minute finite volume corrections to the action, lead to 
errors that are negligible compared to the relevant statistical errors. 

3.2 weak-coupling fits at all N 

We begin by fitting the values of the SU(8) string tension that have been calculated on the 
weak coupling branch in Fig. [2l We attempt to fit all the way to the point at which the 
system tunnels through to the strong coupling phase. Although at lower f3 part of this branch 
is metastable, the barrier to the strong couphng phase is so large that its metastability is 
irrelevant. One can see this from the fact that the distribution of values of the plaquette 
(averaged over the whole volume for a given gauge field) shows no sign of a tail developing 
towards the very different values that it would take in the strong coupling phase. 

Working in the mean-field improved coupling scheme s = /, we find that we can obtain an 
excellent fit to all the values in Tabled using the O(a^), 3-loop truncation of eqn(l5]) 

a^a) = 4f (1 + c'yo) (I + J^) ^ A''. (16) 

The best fit is illustrated in Fig. (Hand the fitted values of \fo j t^i and the constant c/ are 
given in Table [T] We see that cj- = 1.18 ±0.04 is indeed 0(1) as naively expected. In fact any 
attempt to fit without a lattice correction, i.e. with c/ = 0, fails very badly, even if we include 
in our fit only the values of a^a at the weakest couplings. Thus for A^ = 8 it is clear that the 
lattice spacing corrections are important and once taken into account they allow a very good 
fit with the usual 3-loop perturbative running coupling. This turns out to be equally true for 
SU(6), and we list the fitted parameters in Table [T] 

Having established at these larger A^ that lattice spacing corrections are indeed needed, we 
proceed to lower A^ using the same functional form as in eqn (fT6l) . For SU(4), where the bulk 
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transition has weakened to a relatively sharp cross-over [10], a good fit with eqn fll6l) is still 
possible but only if we restrict ourselves to somewhat weaker couplings. This is displayed in 
Table [7] where we express the fitted range in terms of Oa/ct so as to provide a common physical 
measure of comparison for different A^. In SU(3) the bulk cross-over becomes smoother and we 
have to move even further into weak coupling to obtain an acceptable fit. This is illustrated 
in Fig. [5l Finally for SU(2) the bulk crossover is smoother still and we find that even if we 
include only the values at weakest coupling, we obtain an uncomfortably large per degree of 
freedom, as shown in Table [71 All this accords with the naive expectation that the smoother 
the strong-to-weak coupling crossover, the further into weak-coupling one has to go before 
weak-coupling expansions become applicable. 

The computational cost of SU(A^) lattice calculations grows roughly oc A^^, and so our 
SU(3) calculations extend to much smaller values of a than for SU(8). One might wonder if 
it is this greater range that makes the fit to SU(3) so much more difficult than for SU(8). In 
fact this is not the case. If we fit the SU(3) values over a range of a^a values that is very 
similar to that of the whole SU(8) range, we obtain a statistically unacceptable /udf ~ 8. 

We can now take each value of ^Ja j Ki in Table [7] and transform the ratio to \fa / Kjjg using 
eqn llTT]) . If we then extrapolate to A^ = oo using the expected leading 0{1/N'^) correction, 
we obtain a good fit for A^ > 3, giving 

^ = 0.503(2) + ; A^>3. (17) 



The final result for Ajfg is remarkably precise, but this appearance is quite misleading since 
the quoted error is only statistical. The systematic errors, indicated by the '?' in eqn llTTl) are 
potentially very much larger. The most significant uncertainty, the choice of coupling scheme, 
is what we address first. There are some further smaller systematic errors (which are still 
much larger than the statistical errors) which will be discussed later on in this Section. 

One might have hoped hope that the accuracy of our calculation of ay/a would prove 
sufficient to usefully constrain what are the best coupling schemes to use. Unfortunately this 
turns out not to be the case. To show this most graphically we consider three alternative 
schemes which all have the same scale = A/. These are obtained by replacing the value of 
the average plaquette Up in eqn® by a truncation of its perturbative expansion in eqn([9]) to 
either 0{g'^), or 0{g^), or 0{g^). We label these schemes by s = Ji, I2, 13 respectively; i.e. 

1 1 

7^ = 72 (1 - ^1^1 ^J^L ) (18) 

using the notation of eqnQ. For j > 1 these schemes are identical to the mean-field scheme, 
s = /, to 0{g'^), so they will have exactly the same Ag parameter. We now take eqn fll6l) and 
fit our SU(8) string tensions, replacing the values of by the corrsponding values of g\ for 
i = 1, 2, 3. We obtain perfectly good fits in all CcLSGS cLS shown in Table [HI However, as we 
see there, the fitted values of a/ct/A/ vary by almost a factor of 2 between these four coupling 
schemes. (For completeness we show the fit using the straightforward lattice coupling scheme 
s = L, which, when translated to a value of ^/a/Aj^g, is similar to the fit in the I2 scheme.) 
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We obtain very similar results for = 6. Thus in the case of larger N, where the presence 
of a first order bulk transition defines for us where a weak coupling expansion should begin 
to be applicable, there is in practice no useful constraint on the coupling scheme. At smaller 
A^, where there are in fact marked differences in where one can begin fitting with different 
schemes, the presence of a smooth strong-to-weak crossover means that it is not clear what 
significance one should read into any particular comparison. 

3.3 SU(3) : choosing a scheme 

The simplest way to determine if any given bare coupling scheme is 'good' is to compare it 
against a running coupling that has been calculated on the lattice with enough statistical and 
systematic precision to serve as a benchmark. A good example is the Schrodinger Functional 
(SF) scheme of the Alpha Collaboration [27]. In [13] the SU(3) calculation [12] has been 
extended so that it covers a range of energy scales comparable to that covered by experimental 
measurements, but with much greater precision. (Compare Fig.4 of [13] with Fig. 10 of [11].) 
Its (continuum) running is well described by its 3-loop /9-function [T2l [13] . while the size of 
the non-universal P2 term is known to be modest [28l [13]. All this encourages us to use it as 
our benchmark good coupling scheme. 

The SF coupling is a coupling defined on a length scale la for a given a, where / can take 
any integer value. (Restrictions arise in practice.) The values of (3 and / at which calculations 
exist, with the corresponding values of ggp{al), are listed in Tables 3 and 6 of [13]. Our 
strategy is to take some physical quantity /x that has been calculated over a large range of 
lattice couplings, interpolate from the values of (3 at which it has been calculated to the 
values of (3 at which ggp{al) has been calculated for some /, and then fit these values using a 
3-loop perturbative expression, modified by the expected lattice corrections, to obtain fi/AsF- 
(Note that in contrast to extrapolation, interpolation does not generate significant systematic 
errors.) We then fit these same values of /i using eqn( fT6l) with various improved lattice coupling 
schemes to obtain alternative estimates of fi/Ag and hence of h/Asf- If the latter is consistent 
with the value obtained using ggp{al), we have evidence that the corresponding 'improvement' 
is indeed a 'good' one. 

Our string tension calculations in Table [21 do not have enough overlap in f3 with the q'sfQ'^) 
calculations to be useful for this purpose. (Within the relevant range /3 < 6.515, values of 
g^gp have been calculated at only two values of [3.) Instead we turn to the Tq calculations in 
Table [S] which do have a useful overlap. We fit using 



Here there are two lattice spacing corrections. The usual O(a^) term arises from corrections 
to ro(a) etc. while the 0(1/Z^) term arises from lattice corrections to g'spiJ-^) scale 
I X a. As indicated in [T21 [13] one expects the leading correction for the SF scheme to have 




(19) 
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p = 1. However the fact that the couphng is Symanzik-improved means that the dominant 
correction in our range of a might in fact have p = 2 [121 IS] • We shall use fits with both 
values of p taking the difference as part of our estimate of the systematic error. We expect 
(as usual) that Cr ~ 0{1) but we anticipate that dr will be small given that very small scaling 
violations are seen in the step-scaling function in [I2l [T3] .) 

To obtain as strongly constrained a fit as possible, we want to maximise the number 
of values of the running coupling that we fit in the range f3 G [5.70,6.92] where we have 
calculations of a/vQ. We therefore use not only the values of gspila) and ggp{2la) in Table 3 
of [13] but also the values of gsri^^) listed in Table 6 therein. (We exclude the last 4 rows of 
Table 6 since they use a different improvement.) We list in Table |9] the results of our various 
fits. We show separately the results of using different interpolations in obtaining a/rO. As 
expected we find that this makes a negligible difference. We show fits for powers p = 1 and 
p = 2 and see that this makes a small ~ 2% difference. Finally we give an example of a fit 
that uses values in the additional range (3 G [6.92, 7.26] where we obtain the values of a/ro by 
the less reliable process of extrapolation. Again there is only a small change in the fit. 

We note that the region of couplings P G [6.22, 6.91] that we actually use in our fit, 
corresponds to the range of scales 

^ - ^ - ^ X — G [23.9,58.6]. (20) 



AsF aAsF roAsF a 

If we look at Fig. 4 in [13] we see that in this range, the 3-loop formula we use already appears 
to provide a very good approximation to the continuum /^-function. The lattice spacing 
corrections are less well determined but by taking the range of results spanned by both the 
p = 1 and the p = 2 fits, our estimate of this systematic error should be credible. Putting all 
this together we obtain 

3.2(1) — ^ roAjYs = 0.640(20) (21) 



tqAsf 

using AsF — 0.48811A]g^ [12]. If we now compare this with the values in Table [TO] that have 
been obtained by fitting Vq/q using eqn (fT6|) with various bare coupling schemes, we see that 
the Mean-Field scheme produces values that are consistent. It is therefore plausible to adopt 
the latter as our 'good' scheme, while incorporating a systematic error ~ 6% based on the 
difference between the value in Table [TD] and the value ~ 0.66 that one gets at one standard 
deviation in eqn fl2Tl) . 

As a by-product of these calculations we also obtain an updated value for the relationship 
between the two standard scales tq and ^cr: 

r-o^a = 1.160(6)(6) : SU{3). (22) 

The first error is statistical and the second is systematic. It includes a ~ 0.5% error on 
the value of a and a ~ 0.2% error from O(a^) corrections. (See Section 13.51 ) We cannot 
estimate, and therefore neglect, any systematic error on rg. Note that there are no perturbative 
uncertainties here, as we see from the similarity of the various estimates of ro^/a in Table [TOl 
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What we are doing is equivalent to interpolating a/r^ and a^a to common values of /3, taking 
the ratio, and performing a continuum extrapolation with conventional lattice corrections. 
This process is insensitive to the precise interpolation as long as it is smooth. 



3.4 comparing schemes directly 

Here we introduce a method for comparing different coupling schemes directly, without the 
use of any physical quantity such as ro/a or ay/ a. This has the great advantage that it allows 
us to perform comparisons much deeper into weak coupling. 
For a scheme s define the 3-loop perturbative factor 

(23) 

Now we expect for the SF scheme 

laksF^{l + ''j;]Fi^[gl,{al)] (24) 

and for a 'good' improved scheme /' based on the bare lattice coupling 

aKr^{l + da']Fl'[gl{a)l (25) 

up to the various higher order corrections. If we now replace the a? on the RHS of eqn (l25l) by 
the expression for a in eqn fl24p . and if we then take the ratio of the two equations, we obtain 

A/' ° / Fi[gl{a)] [l + c,^ {l + {Fi^[gl,{al)]Y] 

where C2 = c'/A'^p. We can now perform a fit for the constants Cq, Ci and C2 over ranges of jS 
further and further into weak coupling, and see how rapidly Cq approaches the known value of 
Asf/Aj'. The more rapidly it does so, the 'better' we may judge the coupling scheme to be. 

Before proceeding, some remarks. We do not expect exact agreement between the fitted 
and known values of Asf/Aj/ because we are missing the 4- loop and higher contributions to 
the /3-function. Indeed it is precisely the discrepancy that will tell us how 'good' is our I' 
scheme. In addition, over the wider range of g'^{a) values that becomes accessible with this 
method, we may well need to worry about the fact that our supposed constants, ci and C2, 
are in fact power series in g"^. These, and other issues that arise in such a calculation, we shall 
ignore, with the caution that one should therefore regard our calculation as being less than 
fully quantitative. On the other hand, we note that the O(a^) corrections to eqns fl24ll25l) are 
in fact negligible in such a comparison except at the very smallest values of (3, thus largely 
eliminating one source of systematic error. 

We have grouped the calculated values of ggp by the ranges of /5 within which they have 
been calculated, and we have fitted each such group of values separately. Each group has ~ 10 
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values of g^p which is more than enough to constrain the three parameters in the fit. We can 
thus see how the fitted value of Cq compares to the ratio of the A parameters as we go further 
into weak coupling. We do so separately for the mean-field coupling, /' = J, for the variation 
/' = Is where the plaquette is replaced by its 3-loop approximation, and finally for the lattice 
bare coupling, I' = L. We choose to take a power p = 2 although in practice we get much the 
same picture with p = 1. 

The results of our fits are shown in Fig. [61 We plot there the ratio of Cq to its asymptotic 
value, i.e. Asf/-^i for schemes I and J3, and Ksf/^l for the lattice scheme. We see quite 
clearly that Parisi's original mean-field scheme, using the full plaquette in the improvement, is 
consistent with converging remarkably quickly to the expected value (with a small deviation 
that is consistent with a slowly decreasing 4-loop term). If the SF coupling is 'good', then 
so, it would appear, is the mean-field scheme. The other schemes considered fare much less 
well by this criterion. 

A final aside. We have emphasised that the SF coupling possesses various desirable prop- 
erties - it has been calculated over a very large range of scales, it is very precise, it has modest 
higher order perturbative corrections (as indicated by the results of [121 113] )• All this makes 
it an attractive benchmark coupling. However we also recall that Qspi^) is defined on a 
torus. Now, it is known that SU(A^ = 00) gauge theories suffer a sequence of phase transitions 
as / ^ (see [31] and references therein). Each of these transitions involves the breaking of 
a Zn symmetry in one of the space-time directions. The first transition is the usual decon- 
fining transition and the remaining ones can be interpreted as continuations of deconfining 
transitions on ever more dimensionally reduced space-time volumes [32]. At finite N these 
phase transitions will become cross-overs. Such cross-overs will, in general, contribute non- 
perturbatively to the running of g'sp{l)- Although one might argue that such contributions 
will be negligible at the level of accuracy we are aiming for here, this needs to be checked and, 
in any case, this raises interesting issues (in its own right) which need to be addressed. 

3.5 other systematic errors 

There are a number of other, relatively straightforward, systematic errors that we now address. 
These include errors from neglecting higher order corrections in 1// in eqn fllSj) when extracting 
a?a from the loop masses, and also errors in the actual calculation of these loop masses. There 
is the error due to the neglect of higher order corrections in in eqn f|T6l) . Then there is an 
error in using only a leading 0(1/A^^) correction in eqn ffTTl) . Finally there are the errors arising 
from our ignorance of higher order perturbative corrections in eqn fll6l) . within the mean-field 
improved coupling scheme. Assuming this to be a 'good' scheme, as argued above, these last 
corrections can be plausibly bounded. 

We can obtain an estimate of the magnitude of the correction to eqn (fT5|) from the calcula- 
tion in [26]. One finds that for SU(6), and SU(4), the extra correction can be fitted by a term 
that one can rewrite as — 1.23(21)/a^(jL^. For our lattices a?aL'^ ~ 10, so this corresponds to 
a ~ 10% addition to the bosonic string correction, and hence a ~ 0.5% increase in the final 
estimate of a^J(T. This provides us with an estimate of our systematic error from this source. 

We extract the loop masses by identifying effective mass plateaux in appropriate correlators 
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obtained from a variational calculation. In practice the best variational ground state has an 
overlap of ~ 99% onto the true ground state so any error should be very small. In the 
technically very similar case of -D = 2 + 1 that was analysed in [29] , the shift in the extracted 
value of a^a induced by excited string states was shown to be less than 0.5%. We can take 
this as a bound on the corresponding systematic error in the present calculation. Note that 
this error will decrease the string tension and will therefore partially cancel against the error 
discussed in the previous paragraph. 

The fitted values of ^Ja j Ki turn out to be very robust against the inclusion of additional 
O(a^) corrections in eqn (|T6l) . with shifts that are typically ~ 0.25%. 

The perturbative expression in eqn fll6l) is missing an extra factor ~ (1 + Yl^ng'^^) that 
arises from the unknown higher order perturbative corrections. (Ignoring complications that 
arise from the - at best - asymptotic nature of the expansion.) One might imagine that the 
simplest way to proceed is to try and fit the first couple of terms in this series. This turns 
out not to work. The reason is that over our range of scales, the coupling g'^{a) varies very 
little; thus the correction term is to first approximation just a constant (1 + X] '^rifi'^") ~ 
^1 + ^ dng"^^ which renormalises the fitted value of the overall constant coefficient a/o"/A/. 

For example, if we remove from our fit the 3-loop factor exp{— /5|5'|/2/5o} in eqn f|T6|) we find 
that we still get a perfectly good fit, but with a value for y^/A/ that is increased by a factor 
~ 1.14. Now if we use, say, the value g'jN ~ 5.36 that we get at = 44.35 in Table O we find 
that exp{— /5|5f|/2/9Q} ~ 1.17 which is numerically similar. What this means is that neglecting 
higher order perturbative corrections effectively shifts ^/a / away from its true value. In the 
example above, neglecting the 3-loop contribution induces a systematic error in the evaluation 
of i/o"/A/ that is roughly 15%. If we have a well-behaved coupling scheme, one might hope 
that the error due to the neglect of 4-loop and higher order corrections will be roughly the 
square of this, i.e. ~ 2 — 3%. In fact if we estimate, in the above way, the shift induced by 
the 0{g^) term in eqn(l5]) in the MS scheme, where the necessary 4-loop calculation has been 
performed [SU], we find a shift of this magnitude. Assuming that our mean-field improved 
scheme is also a good one, it then seems safe to bound the systematic error from this source 
by ~ ±4%. 

While our 2-loop expression in eqn f|T6|) is exact, the 3-loop expression is kept only to 
leading order in g'^. This makes perfect sense since at higher orders in g"^ one obtains the 
unknown contributions from higher loops discussed above. Nonetheless we have checked the 
difference it makes to our fits and extrapolations if we do employ an exact 3-loop expression 
(by using an expansion to much higher order in g'^) and what we find is that the shift in 
calculated quantities is only at the ~ 0.5% level for both the SF and mean- field coupling 
schemes. 

Finally, we find that the inclusion of an additional 0(1/A^^) correction in the large- 
extrapolation in eqn fll7l) leads to a shift of no more than ~ 0.5%. 

Note that all the systematic errors we discuss will to a first approximation be independent 
of A^ and should therefore not affect the quality of the large- A^ extrapolation, but will merely 
add extra uncertainties to the fitted values. 

We thus estimate a reasonable bound on the systematic error from these sources to be 
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±5%. Adding (in quadrature) the ~ 6% systematic error coming from the choice of scheme 
(as estimated in Section I3.3P we shall take ±8% as our total systematic error. 



4 Conclusions 

As a by-product of the calculations in Section 13.31 we obtained an updated value for the 
relationship between the two standard scales tq and in the continuum limit : r^^a = 
1.160(6)(6) where the first error is statistical and the second is systematic. 

One of the main purposes of our calculation was to determine the importance of lattice 
corrections in the relationship between g^[a) and the lattice spacing a, where g^[a) is the bare 
lattice coupling, or some improvement thereof. The novelty of our approach is to do so for 
SU(6) and SU(8), where the first-order strong-to-weak coupling transition gives unambiguous 
guidance as to where a weak-coupling fit should be applicable. By contrast, in SU(3) the 
presence of a smooth strong-to-weak coupling cross-over makes it difficult to evaluate the 
apparent success or failure of different weak-coupling fits. 

We found that, at larger A^, the variation of a^a with g^{a) did indeed demand the 
presence of an O(a^) correction to the 3-loop perturbative running, and that its coefficient is 
0(1) X cr as one would naively expect. Indeed, such a fit accurately describes the running of 
g^ia) along the whole weak-coupling branch, including the metastable portion that extends 
well beyond the location of the bulk transition. 

As we go to lower A^, we find that we need to go deeper into weak-coupling before such 
weak coupling fits start working. Indeed for SU(2), it is not clear if there is a significant such 
region below /? = 2.70 (the weakest coupling at which we have a calculation of a}a\ An 
interesting question concerns the functional form in the extended cross-over region where the 
strong coupling expansion gradually transforms into a weak-coupling one. As we discussed 
in Section [21 a good place to begin such an analysis is in the D = W \ case where analytic 
arguments are possible [30]. 




Having a well-motivated lattice correction [9] to the running coupling, invites us to address 
the more ambitious goal of calculating K-^/^a for various A^, so as to interpolate and ex- 
trapolate to all A^, including A^ = oo. However in the range of scales a^^a where calculations 
exist, different perturbative schemes give quite different results. This well-known problem is 
usually addressed by 'improving' the lattice bare coupling. However this turns not to be suf- 
ficient. Even if we use variations on the improvement that leave the A parameter unchanged, 
one can easily have a factor of two variation in the fitted Aj-fg/^/a, as we saw with the oth- 
erwise well-motivated mean-field improved coupling. Unfortunately our large- A^ calculations 
turn out not to be able to discriminate between these different coupling schemes - typically 
they all work almost equally well as fits. To finesse this impasse, we turned to SU(3) and 
compared a number of variations on the mean-field scheme with the Schrodinger functional 
(SF) scheme yjj which has been calculated over a range of energy scales comparable to that 
of experimental determinations of ^^(Q^) and with greater precision. Applied to calculated 
values of Tq this comparison indicated that the original Parisi mean-field scheme works well. 
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Using this scheme for all our values of A^, we obtain the fit to all > 3 shown in Fig. [71 

^P = 0.503(2)(40) + ^^:^^^^ ; N>3 (27) 



where the first error is statistical and the second much larger error is expected to provide 
a bound on the systematic error from all sources. (A fit including the N = 2 value is not 
statistically excluded, but we prefer not to include it given the marginal nature of the SU(2) 
calculation.) Note that the component systematic errors are independent of A^ to a first 
approximation, and we can therefore include that error as a common factor to the best fit in 
eqn(l27]). 

We also presented a way of comparing different coupling schemes that does not depend 
on the explicit calculation of any physical quantities. This means that one can perform the 
comparison much deeper into weak coupling. Here again we found evidence that the mean-field 
improved coupling is a 'good' one in the sense of possessing small higher order corrections. 

Our calculations have, of course, have been limited to a small set of improved lattice 
couplings. There will certainly be other 'good' lattice couplings, and the methods described 
in this paper can help to identify them. We have also focussed on one particular lattice action 
(albeit the one that has been most widely used); however it is straightforward to construct a 
Mean-Field improved coupling for other actions, in the spirit of 0. 
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SU(2) 


P 


r 
L 


ami 


Up 


9 1 7(^8 
Z. i ( Oo 


c 
o 




n 99 

U.OOiZZ 


2.2400 


8 


1.380(16) 


0.58286 


2.2986 


10 


1.240(13) 


0.60180 


2.3715 


12 


0.9071(94) 


0.62272 


2.3726 


12 


0.9088(63) 


0.62302 


2.4265 


16 


0.8470(68) 


0.63632 


2.5115 


20 


0.5728(54) 


0.65421 


2.5500 


20 


0.4363(40) 


0.66137 


2.6000 


24 


0.3778(44) 


0.67001 


2.7000 


32 


0.2688(30) 


0.68557 



Table 1: The mass, ami, of a closed flux loop of length L, and the average plaquette, Up, at 
the indicated values of P in SU(2). 



SU(3) 




L 


ami 


Up 


5.6500 


8 


1.425(25) 


0.53750 


5.6750 


8 


1.249(15) 


0.54366 


5.6925 


8 


1.130(12) 


0.54756 


5.6993 


8 


1.106(11) 


0.54896 


5.7995 


10 


0.8860(79) 


0.56755 


5.8000 


10 


0.8766(86) 


0.56764 


5.8945 


12 


0.7283(67) 


0.58111 


6.0625 


16 


0.5408(46) 


0.60034 


6.2000 


20 


0.4465(32) 


0.61362 


6.3380 


24 


0.3588(29) 


0.62560 


6.5150 


32 


0.2943(35) 


0.63948 



Table 2: As in Table □ but for SU(3). 
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SU(4) 


/? 


L 


anil 


Up 


10 


480 


8 


1.348(15) 


0.5253 


10 


500 


8 


1.234(13) 


0.52920 


10 


550 


8 


0.9878(87) 


0.53732 


10 


590 


10 


1.1296(94) 


0.54259 


10 


635 


10 


0.9587(105) 


0.54764 


10 


637 


10 


0.9541(43) 


0.54789 


10 


700 


10 


0.7789(90) 


0.55408 


10 


789 


12 


0.7885(56) 


0.56169 


10 


870 


12 


0.6456(72) 


0.56793 


11 


085 


16 


0.5663(54) 


0.58239 


11 


400 


20 


0.0044(44) 


0.60044 



Table 3: As in Table □ but for SU(4). 



SU(6) 


/? 


L 


ami 


Up 


24.300 


8 


1.246(14) 


0.51721 


24.350 


8 


1.077(11) 


0.52254 


24.425 


10 


1.210(12) 


0.52810 


24.500 


10 


1.063(12) 


0.53258 


24.515 


10 


1.041(11) 


0.53340 


24.670 


10 


0.8431(89) 


0.54089 


24.845 


12 


0.8540(81) 


0.54816 


25.050 


12 


0.6686(60) 


0.55570 


25.452 


16 


0.6396(55) 


0.56866 



Table 4: As in Tabled but for SU(6). 
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SU(8) 


H 


r 




Up 


43.625 


8 


1 279 f 15) 


51220 


43.70 


8 


1.106(12) 


0.51713 


43.78 


8 


0.9989(94) 


0.52094 


43.85 


8 


0.9139(78) 


0.52374 


44.00 


10 


1.0599(91) 


0.52879 


44.35 


10 


0.8036(77) 


0.53848 


44.85 


12 


0.7222(61) 


0.54980 


45.70 


16 


0.6326(45) 


0.56571 



Table 5: As in Table □ but for SU(8). 



SU(3) ; ro 


/? 


ro/a 


Up 


5.70 


2.922(9) 


0.54939 


5.80 


3.673(5) 


0.56778 


5.95 


4.898(12) 


0.58846 


6.07 


6.033(17) 


0.60158 


6.20 


7.380(26) 


0.61384 


6.40 


9.740(50) 


0.63085 


6.57 


12.18(10) 


0.64365 


6.69 


14.20(12) 


0.65204 


6.81 


16.54(13) 


0.65998 


6.92 


19.13(15) 


0.66685 



Table 6: The value of tq [T3], and the average plaquette, Wp, at the indicated values of /3 in 
SU(3). 



N 










2 


[0.177,0.097] 


4.566(28) 


3.83(26) 


1.8 


3 


[0.261,0.101] 


4.888(17) 


2.08(10) 


1.1 


4 


[0.374,0.153] 


5.005(20) 


1.52(5) 


1.0 


6 


[0.415,0.210] 


5.131(25) 


1.30(5) 


0.34 


8 


[0.420,0.209] 


5.199(21) 


1.18(4) 


0.27 



Table 7: Results of weak-coupling fits using eqn( |T6|) over the indicated ranges of ay/a and for 
various A^. 



20 



SU(8) 


s 


ay/a e 


A/A, 








latt 


[0.420,0.209] 


96.80(47) 


3.71(5) 


0.37 


0.3848(19) 


Up 1-loop 
Up 2-loop 
Up 3-loop 
Up meas 


[0.420,0.209] 
[0.420,0.209] 
[0.420,0.209] 
[0.420,0.209] 


8.980(43) 
6.832(43) 
6.100(28) 
5.199(21) 


3.88(5) 
3.67(5) 
3.50(5) 
1.18(4) 


0.54 
0.38 
0.32 
0.27 


0.2932(14) 
0.3854(25) 
0.4316(20) 
0.5064(21) 



Table 8: Weak-coupling fit using eqn fll6p to the whole weak-coupling branch in SU(8); for 
various coupling schemes s. 



SU(3) : SF scheme 



interp. 


P 


Pe 


l/roAsF 




dr 




latt 


2 


[6.257,6.9079] 


3.245(24) 


1.64(80) 


1.05(26) 


0.52 


Up 2-loop 


2 


[6.257,6.9079] 


3.243(20) 


1.61(64) 


1.10(28) 


0.54 


Up 3-loop 


2 


[6.257,6.9079] 


3.238(23) 


1.78(63) 


1.10(28) 


0.38 


Up meas 


2 


[6.257,6.9079] 


3.212(27) 


2.60(55) 


1.13(28) 


0.25 


Up 2-loop 


1 


[6.257,6.9079] 


3.169(23) 


1.63(44) 


0.34(8) 


0.45 


Up meas 


1 


[6.257,6.9079] 


3.141(40) 


2.21(58) 


0.36(09) 


0.18 


Up meas 


2 


[6.257,7.2611] 


3.215(17) 


2.64(45) 


1.09(23) 


0.28 



Table 9: Weak coupling fits to tq in SU(3) using the SF coupling scheme, as in eqn (fT9|) . 



SU(3) : 


s 


a/ro e 


1/roAs 










latt 


[0.166,0.052] 


53.26(21) 


7.03(20) 


1.26 


0.5409(22) 


1,143(11) 


Up 1-loop 
Up 2-loop 
Up 3-loop 
Up meas 


[0.103,0.052] 
[0.166,0.052] 
[0.204,0.052] 
[0.204,0.052] 


5.883(23) 
4.891(20) 
4.567(16) 
4.215(16) 


10.65(30) 
6.93(20) 
5.87(13) 
3.00(13) 


0.89 
1.28 
0.85 
0.48 


0.4475(18) 
0.5383(22) 
0.5765(20) 
0.6246(24) 


1.166(11) 
1.156(14) 
1.159(7) 
1.160(6) 



Table 10: Weak coupling fits to tq in SU(3) using the indicated coupling schemes, s, in 
appropriate modifications of eqn (fT6|) . 
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Figure 1: The SU(3) string tension versus the inverse lattice couphng, including the region of 
the crossover between strong and weak coupling. 
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Figure 2: The SU(8) string tension versus the inverse lattice couphng, including the region 
of the first order 'bulk' transition between strong and weak coupling. Values o are obtained 
coming from strong coupling, while the values • are obtained coming from weak coupling. 
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Figure 3: The (mean-field improved) bare 't Hooft coupling as a function ot the scale in units 
of the calculated string tension, for = 2 (A), = 3 (o), = 4 (*), = 6 (□), and = 8 

(•)■ 



24 




0.1 0.2 0.3 0.4 0.5 



Figure 4: The 't Hooft coupling, defined from tlie mean-field improved lattice bare coupling, 
eqn([8]), as a function of the scale a in SU(8). Shown is the 3-loop perturbative running 
modified by a O(a^) lattice correction. 
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Figure 5: The 't Hooft coupling, defined from tlie mean-field improved lattice bare coupling, 
eqnQ, as a function of the scale a in SU(3). Shown is the 3-loop perturbative running 
modified by a O(a^) lattice correction. 
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Figure 6: Calculated values of R= co/{Asf/-^i'), where cq comes from the fit in eqn (l26i) . For 
the /' = I, •, I' = J3, o, and the I' = L, x, lattice coupling schemes. 
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Figure 7: Calculated values of ^jjjs/V'^ versus l/N"^ with a linear extrapolation to = cxd 
shown. 
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